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The realization that electron localization in dis- 
ordered systems (Anderson localization) jlj is ul- 
timately a wave phenomenon ]2], ^] has led to the 
suggestion that photons could be similarly local- 
ized by disorder This conjecture attracted 
wide interest because the differences between 
photons and electrons in their interactions, spin 
statistics, and methods of injection and detection 

may open a new realm of optical and microwave 
phenomena, and allow a detailed study of the An- 
derson localization transition undisturbed by the 
Coulomb interaction. To date, claims of three- 
dimensional photon localization have been based 
on observations of the exponential decay of the 
electromagnetic wave ^ ||, ^, fj], || as it propa- 
gates through the disordered medium. But these 
reports have come under close scrutiny because 
of the possibility that the decay observed may be 
due to residual absorption ^ [To| , [TT|| , and be- 
cause absorption itself may suppress localization 
JH . Here we show that the extent of photon local- 
ization can be determined by a different approach 

measurement of the relative size of fluctuations 
of certain transmission quantities. The variance 
of relative fluctuations accurately reflects the ex- 
tent of localization, even in the presence of ab- 
sorption. Using this approach, we demonstrate 
photon localization in both weakly and strongly 
scattering quasi-one-dimensional dielectric sam- 
ples and in periodic metallic wire meshes con- 
taining metallic scatterers, while ruling it out in 
three-dimensional mixtures of aluminum spheres. 

In the absence of inelastic and phase-breaking pro- 
cesses, the ensemble average of the dimensionless conduc- 
tance (g) = (G) /(e 2 /h) is the universal scaling parameter 
[fiJI of the electron localization transition jj]]. Here (...) 
represents the average over an ensemble of random sam- 
ple configurations, G is the electronic conductance, e is 
the electron charge, and h is Planck's constant. The di- 
mensionless conductance g can be defined for classical 
waves as the transmittance, that is, the sum over trans- 
mission coefficients connecting all input modes a and out- 
put modes b, g = S a (,T a b (Ref. 13). In the absence of 
absorption, (g) not only determines the scaling of trans- 
mission quantities, such as T a b and T a — SbT^ that we 
will refer to as the intensity and total transmission, re- 
spectively, but it also determines their full distributions 
H , |l5L In electronically conducting samples or in 
white paints, (g) ^> 1 and Ohm's law holds, (g) — N£/L, 



where N is the number of transverse modes at a given 
frequency, £ is the transport mean free path, and L is 
the sample length. But beyond the localization thresh- 
old, at (g) « 1 (Rcfs. 12,17), the wavefunction or classi- 
cal field is exponentially small at the boundary and (g) 
falls exponentially with L. Localization can be achieved 
in a strongly scattering three-dimensional sample with a 
sufficiently small value of £ (Ref. 2), or even in weakly 
scattering samples in a quasi-one-dimensional geometry 
of fixed N, once L becomes greater than the localiza- 
tion length, £ = N£ (Ref. 17). The latter corresponds 
to a wire in electronics or to a waveguide for microwave 
radiation. 

In the presence of absorption, however, (g) cannot 
serve as a universal localization parameter because both 
the small value of (g) and its exponential scaling may 
simply reflect the effect of absorption. In this case, the 
decrease in (g) would represent a weakening rather than 
a strengthening of localization. Moreover, the distribu- 
tions of the intensity and of the total transmission are 
affected by absorption and cannot be simply related to 
(g). Furthermore, it has been argued Q that (g) is not 
a natural scaling parameter because fluctuations in con- 
ductance are so large for localized waves that only the 
full conductance distribution or a parameter reflecting 
this distribution properly expresses the nature of trans- 
port. We find for microwave radiation that for L < £ the 
full distribution of the intensity and total transmission 
normalized to their ensemble averages, s a (, = T a b/(T a i,) 
and s a — T a /{T a ), respectively, can be well expressed 
in terms of a single parameter, the variance of the nor- 
malized total transmission, var(s a ), even in strongly ab- 
sorbing samples [jl9], 20 . This suggests that var(s a ) 
might serve as a localization parameter. Large values of 
var(s a ) would be expected for sharp spectra with widely 
spaced peaks that occur for localized waves. Further, be- 
cause var(sa) — 2/3(5) f° r L ^ £,,L a , where L a is the 
absorption length, and because the localization thresh- 
old occurs at (g) rs 1 in the absence of absorption, wc 
make the conjecture that localization is achieved when 
var(s a ) > 2/3. This localization condition may be ex- 
pressed in a familiar form by defining a new localization 
parameter g' = 2/3var(s a ), which reduces to (g) in the 
absence of absorption in the limit var(s a ) -C 1. Local- 
ization is then achieved for g' < 1 whether absorption is 
present or not. 

In previous measurements of total transmission, 
var(s a ) was found to increase sublinearly with length 
for diffusive waves in strongly absorbing quasi-one- 
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dimensional dielectric samples This raised the pos- 
sibility that the value of var(s a ) might saturate with 
length, and that absorption might introduce a cut-off 
length for the renormalization of transport. Here we 
show that, though the presence of absorption leads to a 
decrease in var(s a ), this appropriately reflects a lessen- 
ing of localization effects. The threshold for localization 
occurs at g' w 1, and for smaller values, g' falls exponen- 
tially with length. 

We now consider wave transport statistics in a 
quasi-one-dimensional geometry. Because L is much 
greater than the transverse dimensions of a quasi-one- 
dimensional sample, energy injected at any point of the 
input is equally likely to emerge at any point of the out- 
put, and modes are completely mixed by the medium. 
As a result, statistical measurements at any point on the 
output surface of the sample or in any of its transmis- 
sion modes yield identical results, and depend upon the 
sample geometry only through the ratio £/L = (g). For 
N 3> f , scattering is locally three-dimensional and wave 
transport may be given a universal description. 

We will first consider the scaling of var(s a ) and its 
connection to localization in strongly absorbing weakly 
scattering dielectric samples contained in a copper tube. 
The role of absorption will be investigated by compar- 
ing these measurements to an analysis of the data that 
statistically eliminates the influence of absorption. Mea- 
surements are carried out in samples of loosely packed, 
1 .27-cm-diameter polystyrene spheres with a filling frac- 
tion of 0.52 within the frequency range 16.8-17.8 GHz. 
In these samples, <«5 cm (Ref. 21), giving (a5m for 
a tube diameter of 5 cm. The exponential attenuation 
length due to absorption is L a — 0.34 ±0.02 m (Ref. 19), 
and the diffusion extrapolation length, which gives an 
effective sample length for the statistics of transmission 
||, L = L + 2z b , is z b « 6 cm (Ref. 21). For L <C £, L a , 
diffusion theory gives var(s a ) — 2L/3£ (Refs. 14,15). 
This result is shown as the horizontal short-dashed line 
in Fig. |l|. Measurements of fluctuations in the spectrum 
of total transmission in ensembles of polystyrene sam- 
ples give the results shown as the filled circles in Fig. [I]. 
These results indicate that var(s a ) increases sublinearly 
with length up to L = 2 m, which was the largest length 
at which accurate measurements of the total transmission 
could be made. 

Values of var(s a ) can also be obtained from measure- 
ment of the intensity T ab . Our measurements of intensity 
and total transmission confirm the predicted relation be- 
tween the moments of the normalized intensity and total 
transmission pl|, (s" h ) = n\(s™). This allows us to relate 
the variance of the normalized transmission to the vari- 
ance of the normalized intensity, which is more readily 
measured in microwave experiments, 
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FIG. 1: Influence of absorption and localization, separately 
and together, on var(s a ) in random polystyrene samples. A 
semi-logarithmic plot of £var(s a ) /L is presented to illustrate 
the measured scaling of var(s a ) over a large range of L, as 
well as various theoretical predictions. The upper and lower 
short-dashed lines represent the two limits of diffusion theory: 
L <C £, L a and i a < I C respectively. The filled circles 
are obtained from measurements of total transmission, while 
the filled squares are obtained from measurements of intensity. 
The open circles are the results of an analysis that eliminates 
the effect of absorption, as explained in the text. The upper, 
long-dashed curve is a fit of these results to an expression in- 
corporating the first-order localization correction to diffusion 
theory. 



Using this relation, we are able to extend our study 
of statistics in random waveguides to greater lengths. 
We measured transmitted field spectra in an ensemble 
of 2,000 polystyrene samples with the use of a Hewlett- 
Packard 8772C network analyzer. The calculated inten- 
sity spectra yield var(s ab ), which gives the corresponding 
values of var(s a ) using equation (1). Values of var{s a ) 
obtained in this way for L < 2 m agree to within 3% 
with those shown as the filled circles in Fig. |l|. The re- 
sults for L > 2 m are shown in the figure as the filled 
squares. They indicate a more rapid, superlinear increase 
in var(s a ) relative to the data for L < 2 m. 

In these measurements, the effect of developing local- 
ization and absorption are intertwined. In order to ob- 
tain the values of var(s a ) that would be measured in 
the absence of absorption, the field spectra are Fourier- 
transformed to give the response to a narrow gaussian 
pulse in the time domain. To compensate for losses due 
to absorption, the time-dependent field is multiplied by 
exp(i/2r a ), where t is the time delay from the incident 
pulse and l/r a is the absorption rate determined from 
measurements of the field correlation function with fre- 
quency shift (^3|. This new field is transformed back to 
the frequency domain, tntensity spectra and the distri- 
bution and variance of intensity are then computed. The 
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intensity distributions are in excellent agreement with 
calculations for diffusive waves (lj, which are de- 
scribed in terms of a single parameter (g) . The values of 
var(s a ) found in this way are shown as the open circles in 
Fig. [I]. A fit of the leading order localization correction 
@, var(s a ) = 2L/3£ + 4L 2 /15£ 2 , to the data gives the 
upper long-dashed curve in Fig. [j] with £ = 5.51 ± 0.18 m 
and Zb = 5.25 ± 0.31 cm. The results are consistent with 
independent determination of these parameters |2l[] . The 
difference between the open and filled circles represents 
the amount by which var(s a ) is reduced, and hence rep- 
resents the extent to which localization is suppressed by 
absorption. 

For diffusing waves, var(s a ) is predicted to fall from 
2L/3£ for L < L a , to L/2£ for L ^> L a (Refs. 24,25), fol- 
lowing the lower short-dashed curve in Fig. 0. Notwith- 
standing the initial drop of var(s a ) from 2L/3£, our 
measurements rise above this curve as a result of en- 
hanced intensity correlation, as L — > £. At L — 5.2 m, 
var{s a ) = 0.6, which is close to the critical value of 2/3. 

To study the statistics of s a for localized waves, we 
examine fluctuations of intensity in strongly scatter- 
ing quasi-one-dimensional samples of alumina (AI2O3) 
spheres. A typical spectrum of s a b obtained near the 
first Mie resonance of the spheres is shown in Fig. ||a. 
The sharp peaks in s a b result from resonant transmission 
through localized photonic states in the medium. The 
distribution function P(s a b) calculated for an ensemble 
of 5,000 samples is shown in Fig. ||b and compared to 
the Rayleigh distribution. The measured distribution is 
remarkably broad, with var{s a b) = 23.5 and fluctuations 
greater than 300 times the average value. The scaling 
of var(s a ) determined using equation (1) at a number 
of frequencies is shown in Fig. |3| We find var{s a ) in- 
creases exponentially once it becomes of order of unity, 
as expected for a localization parameter. 

The availability of a measurable localization parame- 
ter makes it possible to determine the existence and the 
extent of localization in a variety of samples. This is illus- 
trated in measurements of localization in periodic metal- 
lic wire meshes containing metallic scatterers. John has 
proposed that photon localization could be achieved by 
introducing disorder in a periodic structure possessing a 
photonic bandgap p6| . Measurements of transmission 
in the wire-mesh photonic crystal show a low-frequency 
gap [^7j , which fills in as the scatterer density is increased 
|28f . But such measurements leave open the question of 
whether the radiation is localized. To answer this ques- 
tion, we obtain var(s a ) for two concentrations of alu- 
minum spheres shown in Fig. ||. At a volume fraction 
of aluminum spheres of 0.05, a window of localization is 
found, in which var(s a ) > 2/3. At twice this aluminum 
fraction, the reduced values of var(s a ) indicate that wave 
propagation is diffusive. 

We have also used measurements of var(s a ) to exam- 
ine the claim that localization can be achieved in three- 
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FIG. 2: Statistics of intensity in quasi-one-dimensional alu- 
mina samples. Alumina spheres (diameter d a = 0.95 cm, 
dielectric constant e a — 9.8) embedded in Styrofoam shells 
(ds = 1.9 cm, es = 1.05) at an alumina volume fraction of 
0.068 are contained in an 80-cm-long, 7.3-cm-diameter copper 
tube. A typical spectrum of the normalized intensity s a t (a) 
and the distribution P(s a b) (b) near the first Mie resonance of 
the alumina spheres are presented. The sharp and narrow line 
spectra and giant fluctuations shown have been predicted for 
localized waves, and are unlike corresponding spectra and dis- 
tributions in diffusive samples j23j. The distribution P(s a b) 
plotted on a semi-logarithmic scale is computed in an ensem- 
ble of 5,000 sample configurations within the frequency range 
9.88-10.24 GHz, in which statistical parameters do not change 
substantially. The broken line shows the Rayleigh distribu- 
tion. 



dimensional samples of metal spheres at various concen- 
trations |, ||, |oj. We find that in samples of 0.47- 
cm-diameter aluminum spheres of length L = 8.2 cm and 
diameter d — 7.5 cm, with various volume fraction from 
0.1 to 0.475, var(s a ) never rises above the localization 
threshold of 2/3. A maximum value of 0.29 is reached at 
a concentration of 0.45. Thus three-dimensional localiza- 
tion is not achieved in these aluminum samples. These 
results demonstrate that var(s a ) is a powerful guide in 
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FIG. 3: Scaling of var(s a ) in alumina samples. The values 
of var(sa) averaged over the indicated frequency intervals are 
obtained using equation (1). Above a value of the order of 
unity, var(s a ) increases exponentially. In the interval 9.88- 
10.24 GHz, var(s a ) fa exp(L/L exp ), with L exp w 42 cm. 
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FIG. 4: Var(s a ) versus frequency in a wire-mesh photonic 
crystal containing metal scatterers. The photonic crystal is 
a simple cubic lattice made up of copper wires, with a lat- 
tice constant of 1 cm. The lattice has 8 unit cells along each 
side. It is enclosed in a section of a square waveguide and 
filled with mixtures of 0.47-cm-diameter aluminum and Teflon 
spheres, the latter used to dilute the aluminum scatterers. 
Measurements of intensity are carried out in 200 sample con- 
figurations, and var(s a ) is obtained using equation (1). The 
broken vertical line indicates the position of the band edge 
in a periodic structure filled only with Teflon spheres. At 
an aluminum sphere volume fraction (/ a] ) of 0.05, var(s a ) is 
markedly higher near the edge, rising above the localization 
threshold of 2/3 shown as the broken horizontal line. 



the search for and characterization of photon localization. 
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